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Abstract: In 2+1 dimensions, Dirac fermions in reducible, i.e. four-component repre¬ 
sentations of the spinor algebra form the basis of many interesting model field theories 
and effective descriptions of condensed matter phenomena. This paper explores lattice 
formulations which preserve the global U(2 Nf) symmetry present in the massless limit, 
and its breakdown to XJ(Nf)®U(Nf) implemented by three independent and parity- 
invariant fermion mass terms. I set out generalisations of the Ginsparg-Wilson relation, 
leading to a formulation of an overlap operator, and explore the remnants of the global 
symmetries which depart from the continuum form by terms of order of the lattice 
spacing. I also define a domain wall formulation in 2+1+Id, and present numerical evi¬ 
dence, in the form of bilinear condensate and meson correlator calculations in quenched 
non-compact QED using reformulations of all three mass terms, to show that XJ(2Nf) 
symmetry is recovered in the limit that the domain-wall separation L s —> oo. The pos¬ 
sibility that overlap and domain wall formulations of reducible fermions may coincide 
only in the continuum limit is discussed. 
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1 Introduction 

Relativistic fermions moving in two spatial dimensions are at the heart of many interest¬ 
ing issues in theoretical physics, and are increasingly important in effective descriptions 
of phenomena in condensed matter systems. Let us immediately draw a distinction be¬ 
tween so-called irreducible formulations in which the fermion fields if, are described by 
two-component spinors, and reducible models which invoke four-component spinors. In 
the former case, there is no mass term invariant under a discrete parity inversion; gauge 
theories of irreducible fermions generically manifest a parity anomaly mm®, leading to 
an induced Chern-Simons term endowing the gauge degrees of freedom with mass. Such 
theories describe excitations with fractional statistics, and form the basis for effective 
descriptions of the fractional quantum Hall effect [4], The main focus of this paper, 
by contrast, will be reducible theories, which like their 4 d counterparts admit parity- 
invariant fermion masses, but whose action is invariant under an enlarged group of global 
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symmetries generated by both y 5 and the “unused” 73 . The prototype model involv¬ 
ing reducible fermions is QED 3 , an asymptotically-free theory with potential non-trivial 
critical dynamics in the infra-red, and still investigated as a toy model of walking tech¬ 
nicolor [5]. Reducible fermions figure in effective descriptions of the spin-liquid phase of 
quantum antiferromagnets M, the pseudogap phase of cuprate superconductors M, 
and graphene unumna- 

Once interactions are introduced, reducible 3d fermions may exhibit spontaneous 
mass generation via formation of a bilinear condensate (V^) 7 ^ 0 , entirely analogous to 
chiral symmetry breaking in QCD. This is manifested in variants of the Gross-Neveu 
(GN) model [13], in which the basic interaction is a four-fcrmi contact between scalar or 
pseudoscalar densities; and also in the Thirring model [14] whose interaction is a con¬ 
tact between two conserved currents, and which shares the global symmetries of QED 3 . 
For GN models the condensate may be studied systematically using an expansion in 
1/Nf where Nf is the number of fermion flavors [15], whereas condensate formation 
in the Thirring model is inherently non-perturbative. In both cases the critical cou¬ 
pling at which the condensate forms is associated with a UV-stable fixed point of the 
renormalisation group, although the nature of the fixed-point theories remains an open 
problem |l 6 j . 

Because field fluctuations near a fixed point can be large, it is natural to explore 
critical dynamics in 3d using lattice field theory techniques. We immediately confront the 
issue of how best to formulate the fermion fields. The two traditional approaches which 
mitigate the species doubling present in naive discretisations are Wilson fermions, which 
eliminate unwanted species by explicitly breaking important global symmetries, and 
staggered fermions, which preserve a subgroup of the symmetries by effectively spreading 
the spinor degrees of freedom over several lattice sites (for an excellent introduction to 
lattice fermions see chapters 5,7 and 10 of m- In a seminal work [18] Coste and Liischer 
argued that the Wilson formulation is the more natural for irreducible descriptions, 
and can be shown to recover the correct parity anomaly (a recent calculation of the 
mass-dependence of the anomaly in various background fields using Wilson fermions is 
given in [H]), while the staggered approach naturally leads to a parity-invariant mass 
term; indeed the relation between staggered holds and the Nf = 2 reducible spinor 
flavors expected in the continuum limit was given by Burden and Burkitt in [20j. Since 
the Chern-Simons action is imaginary in Euclidean metric, numerical simulations of 3d 
fermions have to date been restricted almost exclusively to this second case. 

To date there have been several papers studying critical behaviour in both GN 
[m eh Eg and Thirring [23[ 123, 25] models using numerical simulations of staggered 
fermions. In general the results support the theoretical prejudice that GN models ex¬ 
hibit a critical point at a coupling g 2 ~ 0(a) (where a is lattice spacing) for all Nf, 
with small corrections of 0(1 /Nf) to both g 2 and the “mean held” critical exponents 
predicted in the large- Nf limit [T5] , By contrast, the Thirring model exhibits gap for¬ 
mation only for Nf < Nf c , with both exponents and g 2 strongly dependent on Nf [ 25] . 
In the continum, the two models are supposedly distinct. However, recent results ob¬ 
tained with a fermion bag algorithm permitting simulation in the massless limit [20, [27]] 
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suggest that in fact the lattice models defined using staggered fermions may actually lie 
in the same universality class for the minimal case Nf = 2. This somewhat unexpected 
result is motivation to question the applicability of staggered fermions to problems in¬ 
volving critical fields. For N staggered flavors the global symmetry is XJ(N)®XJ(N), with 
Nf = 2 N, broken to U (N) by the generation of a fermion mass. In the long wavelength 
limit ad —* 0 this enlarges to the required U(4iV) —>-U(21V)<g)U(2iV) [2D], but there is 
less reason than usual to trust this restoration once fluctuations on all length scales are 
important. Other situations where it may be important to reproduce the global symme¬ 
try pattern correctly include the infrared behaviour of QED3, where it has been argued 
that dynamical mass generation depends on the correct counting of massless degrees of 
freedom, which include Goldstone modes in a symmetry-broken phase [28], and the role 
of “half-instantons” in the thermal response of the gapped phase of graphene, which 
may result in a much reduced value for the critical temperature [2D]. 

This paper explores the application of formulations originally developed to optimise 
the reproduction of global symmetries in lattice QCD, namely Ginsparg-Wilson (GW) 
fermions [30] and, principally, domain wall fermions [31, 32], to reducible fermion mod¬ 
els in 2+1 d. After reviewing the relevant symmetries and identifying three distinct but 
physically equivalent formulations of the mass term in the next section, in Sec. [D] we 
generalise the GW relation to fermions in 2+1 d and identify remnant quasi-global sym¬ 
metries, which recover the desired U(2iV/) form only in the continuum limit a —> 0. 
A realisation of the GW symmetries by an overlap operator [33] is given. In Sec. [4] 
we define a domain wall fermion operator in 2+1+ld which permits the definition of 
fermi fields localised on domain walls at either end of the newly introduced 3 direction 
which purport to satisfy the XJ(2Nf) symmetry in the limit that the wall separation 
L s —> 00 . An important component of the argument is the reformulation of the three 
distinct mass terms given in Sec. [2] Sec. [5] presents results from numerical investigations 
of the Nf = 1 domain wall operator in the context of quenched non-compact QED 3 , 
which permits the use of either weak, strong, or intermediate coupling. While there 
is no attempt to explore either continuum or thermodynamic limits, we calculate both 
bilinear condensates (Sec. 15.111 and meson propagators (Sec. 15.21) using each of the three 
alternative mass terms, and show that in almost all cases as L s —> 00 the results are in 
accord with a scenario in which U(2) symmetry is broken to U(l)<g)U(l). Interestingly, 
the most rapid convergence to the U(2)-symmetric limit is obtained for the case of a 
“twisted” mass term For intermediate coupling the results for the conden¬ 

sate ( , 0'0) are compatible in the massless limit with old results obtained with staggered 
fermions [35]. Finally in Sec. [6] we present a summary of the findings and an outlook 
for future investigations. We also discuss the intriguing possibility that for reducible 
theories of fermions in 2+1 d the overlap and domain wall approaches may not coincide 
except in the continuum limit. 
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2 Relativistic Fermions in 2+ld 


I begin by reviewing the continuum formulation of a gauge theory with fermion fields 
T, 0 in a reducible representation of the spinor algebra, based on 4 x 4 Euclidean Dirac 
matrices with { 7 ^, 7 *,} = 26^, n, v — 0 , 1 , 2 , and having a parity-invariant mass. 
The weakly-interacting long-wavelength limit of staggered lattice fermions naturally re¬ 
produces this formulation with Nf = 2 flavors [20] - in what follows flavor indices are 
suppressed. The action can be written (for convenience, the necessary J d 3 x is omitted 
in all action definitions) 

S = TDT + mTT ( 1 ) 

where the covariant derivative operator D can be expanded as 

D = 7 0 T>o + 7i£>i + 72 D 2 = -D t . (2) 

This has global symmetries 

T H- e ia T ; T ^e~ ia , (3) 

Vi; e *7375 vI ; . ^ ^ ^r e - a7375 , (4) 

where 73 and 75 are two additional traceless, hermitian, and linearly independent 4x4 
matrices which anticommute with all the 7 ^ (see (1K1I91) below), and as usual in Euclidean 
rnatric 75 = 70717273 - For fermion mass m — 0 there are two additional symmetries 

$ K e” 75 4 ; $ fe” 75 , (5) 

T i->- e ia 73 \l> ; T i-> Te*" 73 . ( 6 ) 

These four rotations generate a global U(2) invariance, which generalises to U(2iV/) 
for several flavors. The mass term explicitly breaks the symmetry from U(2i\M —> 

ft will prove interesting to explore different forms of the mass term, which are simply 
accessed by changing integration variables in the path integral. Since there is no axial 
anomaly in 2 + Id, this procedure is straightforward in the continuum and the resulting 
action describes identical physics. If, however, the representations of the Dirac matrices 
are tied to the particular form of the underlying lattice, as is the case for staggered 
fermions or graphene, then due to discretisation effects the mass terms are not equivalent 
and correspond to distinct patterns of symmetry breaking (see the discussion following 
eqn. (Eli for an example). Let’s recast the continuum action (JT|) in terms of two two- 
component spinors u and d: 

S = uDu — dDd + muu + mdd. (7) 

where D = = criDo + oy-Di + (T 3 D 2 and the cq are Pauli matrices. The link with 

dU) requires the identification 
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implying 



We now define an important discrete symmetry, parity, here specified for convenience in 
terms of reversal of all three spacetime axes t—>■ —x^ (in general parity must invert 
an odd number of axes, since flipping an even number is equivalent to a rotation: the 
Euclidean parity operation which flips just one axis is formally equivalent to the time- 
reversal operation frequently discussed in condensed matter physics). In fact it can be 
realised in two ways: 

u i — y d] — u] u i — ^ d] d i — ^ — u] i.e. T (->• ry 3 T; 0 H > —IT73 (10) 

u h -y —id] d (->■ — iu ; u^r id] d (->■ iu] i.e. T (->■ 0 (->■ —id/ 75 ( 11 ) 


This should be no surprise, since both 73 and 75 behave identically with respect to the 
7 m appearing in (JTJ) . In either case the parity operation effectively exchanges the u and 
d fields, absorbing the sign change of D under x —> —x, but keeping the mass term 
invariant. 

Now consider a change of basis 

^ = -^(u + d)] X=-^(-w + d); $ = -^(u~d)] y=-L(u + J) ( 12 ) 
in which the action (JTJ) reads 

*=<* *> e 1) (t) = ^ a (f (t )' < 13 > 

The parity transformation leaving (TT3l) invariant is equivalent to (TTOl) . In terms of four- 
component spinors the action (TTH1) can be written 

S = ^(D-im'ya)^. (14) 


Although the mass term imM> r )^> is still parity invariant, it is now an antihcrmitian term 
in the Lagrangian, in contrast to the hermitian mass of (pD) . Superficially it resembles the 
twisted mass sometimes used to formulate lattice QCD in 4 d [36], though the absence 
of an anomaly in 2 + Id permits the term to be flavor singlet. We can also consider a 
different change of basis 


1 i — 1 — i — 

^ = ~^(u + id)] u j = —(u~id)] ip = —(u + id)] u!=-^(u-id) (15) 


yielding 


S=(ij> to) 


D —im\ f 0 
—im —D ) \ uj 


= w) 


D —im\ f -0 


—im 


Dt 


Ld 


(16) 
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with this time a parity transformation fjlip . In terms of four-component spinors (fIBD 
can be written 

S = ^(D (17) 

again, the mass term is parity invariant and antihermitian. 

While the three actions (1X1114117[1 must be equivalent in the continuum, when derived 
from a lattice system such as a tight-binding model of graphene the terms correspond to 
physically distinct gapping instabilities at the Dirac cones [57] . The mass term mSh of 
(IT]) corresponds to a charge density wave in which electrons preferentially sit on one of 
the two sub-lattices A or B [38], whereas a linear combination of m^Ss (fill) and 
(TT71) yields a bond density wave in which electrons are distributed on both A and B 
sublattices in a Kekule texture [39] . 

Note that a mass term proportional to the bilinear 473754 , whether hermitian or 
antihermitian, is qualitatively different. I 11 terms of two-component spinors the resulting 
action reads 

s^i)( D+m («) 

since all elements are proportional to the combination D + m there is no parity operation 
realisable as a linear transformation on 4,4 which flips the sign of D but leaves the 
mass term invariant. This corresponds to the non-time-reversal invariant “Haldane 
mass”, realised in graphene-like systems by alternately circulating currents in adjacent 
half-unit cells m- 

In summary, there are three linearly independent, physically indistinguishable parity- 
invariant mass terms available for continuum four-component Dirac fermions in 2+1 d. 
We will see that this furnishes a non-trivial test for lattice fermion formulations, such 
as the domain wall formulation presented in Sec. |4] in which the matrices 73 and 75 
appear in inequivalent ways. Of course, the transformations (j 12 Ml 5 f) are merely special 
cases of the rotations (I6ll5l) with a — |, so the physical equivalence of the mass terms is 
a consequence of U( 2 ) symmetry. 

3 Ginsparg-Wilson Relations and the Overlap oper¬ 
ator 

The Nielsen-Ninomiya theorem mi m\ famously asserts the impossibility of a lattice 
fermion formulation with a physical (ie. undoubled) spectrum which simultaneously 
respects locality, unitarity (or reflection positivity in Euclidean metric) and the exis¬ 
tence of a conserved axial charge, expressed in four dimensions via the anticommutator 
{D, 75 } = 0. Ginsparg and Wilson [30] proposed a minimal modification to these criteria 
for lattice chiral fermions to be viable, namely that chiral symmetry now be constrained 
by the GW relation 

75 D + D 75 = aD^D, (19) 
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where a is the lattice spacing. The non-vanishing right hand side of (1191) is effectively an 
0 (a) contact term for the anticommutator of the propagator D _1 with 75 , which should 
become unimportant in the long-wavelength limit aD —> 0. 

The GW relation appropriate to odd-dimensional fermions in irreducible represen¬ 
tations of the spinor algebra was investigated by Bietenholz and Nishimura [43], who 
found that the correct parity anomaly is recovered due to non-invariance of the fermion 
measure under a “generalised parity” transformation. For the reducible 2+1 d theories 
of interest in this paper the GW relations generalise to: 

7 5 D + T>7 5 = aD^D] 

7 3 D + L>73 = aD~/ 3 D] (20) 

7375/1 - -D7375 = 0. ( 21 ) 

In addition we require the hermiticity relations 

73-D73 = 75-D75 = D ] . (22) 

The GW symmetries (119H21 j) are realised by the overlap operator [33] 

D ov = a -1 (l + 7 3 sign[ 7 3 A]) = a _1 (l + 7 5 sign[ 7 5 A|) = a -1 ( 1 + f , ) , (23) 

where the matrix A = aD Wds — 1, where D Wlls is the lattice Dirac operator for ortho¬ 
dox Wilson fermions. It has the properties 75 A 75 = 73 A 73 = and * 7375 / 1*7375 = A, 
implying D ov + D°A = aD ov D°A = aD°AD ov . It is manifest that 73 and 75 are treated 
identically in the overlap formalism. A distinct overlap operator, appropriate for irre¬ 
ducible fermions and reproducing the correct parity anomaly, was considered in pH] . 

Analogously to the situation in Ad BSj the GW relations (I19H21|) admit the following 
remnant symmetries, which recover the desired XJ(2Nf) symmetries ([3j[6]) in the contin¬ 
uum and long-wavelength limits a —» 0, aD —)■ 0: 

e ( ia ')' 5 (i- s r))^ . vp (_)• #e(* a ( 1_2 ^) 75 ) 

vp h-> e (* Q73 ( 1 - £ ^))+ : T t-7 vfeM 1 - 2 ^)^) 

T h+ e _ " 7375 T ; T ^ ^e“ 7375 . (24) 

For theories such as the GN and Thirring models defined near a critical point it is not 
a priori clear that continuum and long-wavelength limits coincide. 

To construct mass terms, we define projection operators of two kinds 


P 5± = “(I ±75); 

P 3± = gl 1 ± 73), 

(25) 

1 , 

1 , 

(26) 

P 5± = 2 ^ ± 75); 

p 3± = ± 73 ), 
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with 


75 = 75(1 -glD); 73 = 73(l-a-D)- (27) 

Relations (1191121(1 may be used to show 75 = 7 3 = 1, so that the P s are genuine projec¬ 
tors. Importantly, it also follows that 

DP 3 ± = P 3t D- DP 5 ± = P§zfD. (28) 

Therefore it makes sense to define two different sets of projected fields: 

= Pa±V', ^ 3 ± = ^P3P, 

^5± = Ai'f; ^ 5 ± = ^P 5t , (29) 

enabling a decomposition of the kinetic term in two different ways: 

TUT = T 3+J DT 3+ + T 3 _.DT 3 _ = T 5 + _DT 5+ + f 5 _DT 5 _. (30) 

The mass term in this approach is defined in terms of projected fields: 

m h (V 3 _T 3+ + T 3+ T 3 _) = m h (T ;5 -T 5+ + ^ 5 +^ 5 -) = m h V( 1 - ^)T = m h Sg w . (31) 

Although this term only corresponds with the desired mass term in the continuum limit 
a —» 0 , the extra 0 (a) piece has the same form as the kinetic term, so the overall effect 
is a benign wavefunction renormalisation. We can also check the effects of the GW 
symmetries (1241) . along with the fermion number symmetry (J3]) . With obvious notation, 
and to 0(a ): 


<5f w T(l - — )T = 6g w V(l-— )T = 0 (32) 

^*(1-^)* = ^(1 - 2™>)* (33) 

The variation (133]) provides the interpolators for the Goldstone modes associated with 
a spontaneous breaking (T(l — ^)T) ^ 0; we see there is an 0(a 2 ) correction to the 
expected continuum forms. 

It is interesting to repeat this exercise for the two other parity-invariant mass terms 
in (11411 171) . The mass terms constructed from the projected fields are 


im 3 (T 3 _ 7 3 T 3+ + T 3 + 7 3 T 3 _) = im 3 (T 5 _ 7 3 T 5 _ + T 5 + 7 3 T 5+ ) 
im 5 (T 3 _ 75 ^ 3 - + T 3 + 7 5 T 3+ ) = im 5 (f 5 _ 7 5 T 5+ + T 5 + 7 5 T 5 _) 
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7m 3 T7 3 (l 

m 3 Sg w - 

75(1 

m 5 Sg w . 


aD 

~T 


)T 


aD 

~T 


)* 

(34) 




In this case the 0(a) correction does not correspond to an existing term in the action, 
so the correction may be less benign. Another symptom is that if we label the full Dirac 
operator including mass terms as V = D + rrihSi + + ITI 5 S 5 , then 

'y 3 V(m hl 0, 0)73 = V ] (m/,.,0,0); ^V{m h , 0, 0)75 = T>\m h , 0, 0), (35) 


but that no such simple relations exist once 7773 , 777.5 7 ^ 0 . 

We find the O(a) variations of the mass term S 3 under (1211) 


S? w 473(1 - — )<P 

<S? W *73(1 - 
if W t T3 (l - 
^* 73(1 - 


0 


27aT(l 


a 2 DD' 

4 


)(1 


7a^7 3 75^(D - D ] ){ 1 
- 2 a% 5 (l - ^)tt 


aD 


2 
aD 






(36) 


with similar results, mutatis mutandis for variations of S 5 . In all cases the expected 
expressions for Goldstone interpolators are recovered in the continuum limit. However 
for a > 0 there are subtle differences in the effect of the GW variations on the different 
mass terms; for instance, the two Goldstone interpolators “ 7 / 7 /f’ and resulting 

from 8 ^ w , dyg in ()36|) have realisations differing by 0(a 2 ) away from the continuum 
limit. I 11 this sense, full U(2) symmetry is only recovered by overlap fermions as a —» 0. 


4 Domain Wall Formulation 

In this section we set out an alternative route to undoubled U(27V/)-symmetric fermions 
using the domain wall approach first introduced by Kaplan [3T] and put into the present 
form by Furman and Shamir [32] . Here we follow the treatment set out in [IT]. Define 
4-spinor fields ^(x, s), ^(x, s) living on a 4 d lattice where x denotes the usual 2+ld 
spacetime coordinates of a lattice site and s = 1,..., L s its coordinate along the extra 
dimension, here labelled 3. The kinetic term in the action is then 

S™ = ^2^2^>{x,s)D DW (x,s\y,r)^(y,r). (37) 

x,y s,r 

with domain wall Dirac operator 

D dw (x, s\y, r) = S Str D(x\y ) + S Xjy D^ w (s\r). (38) 

The first term is the orthodox 2 + 1 d Wilson operator 

D { x \v) = \ [( 1 ~^) U ti( x )Sx+ii,y + (l+'Y IJ ,) U l(y)Sx-(i,y]+{ M -3)Sx,y (39) 

z 74=0,1,2 
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with gauge link variables U^x), and D^ w controls hopping in the 3 direction: 

D 3 ^ (sis 7 ) = — [(1 — 73)<Wl,s'(l — &s',L a ) + (1 + 73)^s-l,s'(l ““ ^s',l) — 2(5 s y] . (40) 

Note there are Dirichlet boundary conditions imposed in direction 3, at s — 1 and 
s = L s . The inclusion of D™ explicitly destroys the equivalence of 73 and 75 in the 
dynamics described by the action (157j) . so it will be important to test whether and how 
this is recovered in practice. 

The key idea |3TJ is that the dynamics generated by (fTTl) and (jlUj) . with suitably 
chosen M, results in fermion zeromodes localised on domain walls at s = 1,L S , which 
are also respectively + eigenmodes of 73 . The 2 + 1 d physics we wish to describe is 
formulated entirely using these localised modes (the Wilson terms in (13011401) render 
the would-be zeromodes due to unwanted doubler species non-normalisable in the limit 
L s — y 00 [33]). In particular we need to define 2 + Id fermion mass terms corresponding 
to their continuum counterparts in ([!]), (1T51) and f|17|) . To this end, define fermion fields 
- 0 (x) living in 2+1 d: 


i>(x) = P^(x,l) + P + ^(x, L s )] 

'ip(x) = T(a;, L S )P_ + T(a;, 1)P + , (41) 

where from now on P± = |(1 ± 73). We thus consider actions of the form (1571) supple¬ 
mented by three alternative mass terms: 

m h S h = m h iii/j = m h [$(x,L s )P_'f>(x, 1) + ^(x, l)P+T(x, L s )}- (42) 

m 3 S 3 = im 3 ^ 73 ^ = im 3 [+ (x, L s )^f 3 P_^(x, 1) + f (x, 1 ) 7 3 P+T(x, L s )\; (43) 

m 5 S 5 = im 5 ^ 5 '0 = im 5 [V(x,L a )'y 5 P + '&(x,L 8 ) + T(x, l)^P-^(x, 1)]. (44) 


It is interesting to note that Sh has the same form as the fermion mass term for domain 
wall formulations of 3 +Id physics, and couples fields from opposite walls; S 3 also couples 
opposite walls, but S 5 couples fields living on the same wall. 

I 11 the next section we will examine the numerical consequences of the three terms 
(I42H44I) and in particular check whether they yield compatible, U(2)-symmetric results 
in the L s —» 00 limit. 

5 Numerical Results 

In order to explore the domain wall action (1571) supplemented by one of the mass terms 
(I42H44I) we have performed quenched simulations of non-compact QED, so that U tl {x) = 
exp(i9^(x)) with the real link variable 9 equilibrated using the Maxwell-like action 

SncQED — ^(A+^(x) - A+d M (x)) 2 ; fj,,u = 0 , 1 , 2 . (45) 
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The dimensionless coupling j3 = (g 2 a) _1 , where the dimensionful fermion charge g is 
the natural scale with which to define physical observables. The continuum limit lies at 
/ 3 —> oo. In 2 + 1 d QED is an asymptotically-free theory whose infra-red behaviour is 
still imperfectly understood. Since the non-compact lattice formulation is non-confining, 
numerical simulations are plagued by the slow fall-off of the photon propagator (oc r -2 
for the quenched theory), and the thermodynamic limit is extremely difficult to achieve 
at weak coupling. Quenched simulations with staggered lattice fermions presented in 
[35] suggest that chiral symmetry is spontaneously broken in the ground state in the 
continuum thermodynamic limit, namely (in lattice units a = 1 ) 

lim lim lim ( , 0'0(' m )) 7^ 0- (46) 

m/3—>-0 L//3—>oo (3—>oo 


Lattice sizes up to 80 3 were used in support of this claim. 

Away from the continuum limit massless staggered fermions have a manifest U(l)< 8 ) 
U(l) global symmetry which spontaneously breaks to U(l) either when an explicit mass is 
introduced or a chiral condensate of the form develops. Only in the continuum limit 
can this pattern enlarge to the expected XJ(2Nf) with Nf = 2 [ 20 ]. 

The main aim of this study is not to verify (1461) for domain wall fermions, but rather 
to test the restoration of the correct symmetry breaking pattern, with Nf = 1, as 
L s —> oo. To this end we have performed simulations on system sizes 24 3 x L s , with 
couplings f3 = 0.5, 1.0 and 2.0, corresponding respectively to strong, intermediate and 
weak coupling regimes. Throughout the domain wall fermion mass M in (l39j) was set to 
0.9 (physical quantities are expected to be M- independent; reflection positivity requires 
0 < M < 1), and unless otherwise stated the masses mh,ms,rri 5 in (I42H44|) to 0.01. For 
convenience a hybrid Monte Carlo algorithm was used to generate the {9} ensembles, 
with 100 trajectories of average length 1.0 between measurements; a much more efficient 
Fourier-space method was employed in the original study [35] . 

5.1 Bilinear Condensates 

First we explore bilinear condensates generically defined by (^T^) = d\n Z/dirii, where 
rrii G {nih, m 3 , These form in response either to m; 7 ^ 0 or in the thermodynamic 
limit to strong dynamics. For the latter case in the limit m, —> 0 all three condensates 
should be equal if U(2) symmetry is manifest. It is not obvious that this will occur 
for the action (EZD, both because of the Wilson terms inherent in (I39ll40|h and because 
73 and 75 do not appear in (137|) on an equal footing. By hypothesis, however, U(2) 
symmetry should be recovered as L s — > 00 . 

To begin, we present results obtained using a spatial point source on a conhguation 
generated at f3 = 0.5 (in fact, the numbers result from averaging over 10 spatial sources); 
the systematics are easiest to expose at the strongest coupling. Note from (I421I44[) 
that each condensate gets contributions from two terms: for an< 4 *(^ 73 ^) the t wo 
terms arise from four-dimensional propagators running from s = 1 to L s and L s to 1 
respectively; for ilfip 75 ^) each contribution is from a propagator starting and ending 
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on the same domain wall. Within the working numerical precison each contribution 
is the complex conjugate of the other, so the sum is real. However, it turns out the 
imaginary component parametrises the approach to the U(2)-symmetric limit. Define 
z('l r (l) 73 d'(L s )) = f ('073 '?/’)l s + iAh(L s ) (where the first term is real, and the spatial 
coordinate x is suppressed), and then write: 



% ~ 

2 (0730)-Lg—>-oo T A/ l (L s ) T 6/),(T s ), 

(47) 

(0730 )L a 

X — 

= + e 3 (L s y, 

(48) 

(0750 )l b 

% — 

= ^'fa^Ls^oo + e 5 (T s ). 

(49) 


The residuals A/j and e* must vanish for a U(2)-invariant limit. 
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Figure 1: Residual errors as a function of L s for bilinear condensates evaluated using point spatial 

sources on 24 3 at /3 = 0.5. 
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Figure 2: Residual error 2A h(L s ) evaluated for various /3 using a stochastic estimator on 24 3 . 


Fig. Q] plots the residuals for L s — 16,... ,40; note that A h is measured directly as 
the imaginary component of f(' 073 ' 0 ) using just the + components of \1/, while to 
estimate the e* the value of l s ^oo is taken to be that measured at L s = 48. 

Several features are apparent: 

• The dominant correction by almost an order of magnitude is Ah, which contributes 
to the hermitian condensate but not, as a result of the twist, to the antiher- 
mitian *(^ 73 , 5 ^)- Indeed, at the weakest coupling (3 = 2.0 Ah is the only residual 
large enough to measure. 

• All residuals decay approximately as exp(— cL s ), consistent with U(2) symmetry 
restoration as L s —* 00 . 

• 63 and 65 are practically identical - the disparity in Fig. Q] at L s = 40 probably 
results from uncertainty in the determination of i ('073'0)l s ->oo- This is striking 
since the underlying structures in terms of four-dimensional propagators are quite 
distinct. 

Of course, we should not draw universal conclusions from a single gauge configu¬ 
ration, particularly since the bilinear condensates display strong intermittent upward 
fluctuations across a gauge ensemble. Fig. [2] shows A h(L s ) evaluated for 1200 (24 3 ) or 
2400 (12 3 ) measurements of i(xp 73 ^) at the three couplings explored (for L s > 16 0(500) 
measurements were made at the strongest coupling f3 = 0.5 on 24 3 ) using a stochastic 
noise vector 77 ( 2 ;); to faithfully implement the projectors P± in (14H|) the same 7] is chosen 
for all spin components in the trace. Reassuringly, for L s > 8 the exponential behaviour 
A h oc e~ cLs persists, with decay constant c decreasing as the coupling increases. Com¬ 
parison of the f3 = 0.5 results also shows c is in general also volume-dependent. Although 
not investigated here, we also expect c to depend on the Lagrangian parameter M. There 
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is no reason to doubt, however, that the residuals will eventually become an insignificant 
systematic as L s —» oo. 



Figure 3: Bilinear condensate f° r various /3 as a function of L s . Results for (’ip’ip) and 

at /3 = 1.0 are also shown. 



Figure 4: Average number of conjugate gradient iterations per measurement required for the data of 

Figs. H31 


Results for the condensate which as a result of the twist does not include 

effects of the dominant residual Ah, are shown in Fig. [3l With the use of a stochas¬ 
tic estimator a conjugate gradient residual of 10~ 12 per (four-dimensional) lattice site 
and spin component sufficed. As a function of L s the results plateau essentially once 
Ah/ 1('073'0)| is smaller than the statistical error confirming the trends of Fig. [2j this 
occurs for L s > 8 for /3 = 2.0 and L s > 16 for /3 = 1.0. At the strongest coupling /3 = 0.5 
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the available numerical resources have not permitted this regime to be probed on the 
reference 24 3 volume; however data taken on 12 3 eventually plateau for L s > 28. Based 
on the Ah data from Fig. [2] we might guess the j3 = 0.5, 24 3 plateau may set in for 
L s 40. 

Fig. [3] also shows results for the other bilinear condensates ( , 0'0) (h) and *(^’ 75 ^’) 
(5) for f3 = 1.0, showing that within statistical errors they become consistent with 
*('07 3 '0) for L s > 20. The impact of the A h residual on {'ip'ip) is clearly discernable, and 
suggests it is not the optimal choice for systematic studies of chiral symmetry breaking, 
ffowever, with this level of precision there is no reason not to suppose U(2) symmetry is 
eventually restored. Fig.HJshows the number of conjugate gradient iterations required for 
each fermion matrix inversion is roughly constant once the plateau is reached, implying 
that the numerical effort needed scales linearly with L s as the U(2)-symmetric limit is 
approached. 

Finally, Fig. [5] plots 73 ^} as a function of m 3 at /3 — 1.0 on 24 3 x 20. The data 
can be plausibly extracted linearly to the chiral limit m 3 —y 0 to yield a non-vanishing 
intercept, consistent with the hypothesis that chiral symmetry is indeed spontaneously 
broken (though at this stage no thermodynamic or continuum limit is taken). Also 
shown is the staggered fermion condensate \{xx) (l e - normalised to Nf — 1) on the 
same spatial volume estimated from Fig. 7 of [35]. The values are close enough to 
be plausibly consistent, but in any case provide reassurance first that the domain wall 
fermion data do not contain significant contributions from doubler species, and secondly 
that the explicit chiral symmetry-breaking violation by the Wilson terms in (IHUjl is 
mitigated in the domain wall approach. Of course, since there are as-yet unquantified 
renormalisations of both elementary fields and composite operators, these considerations 
strictly require a continuum limit to be taken. 
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7T 




• L s = 20 
0-0 staggered 


0.01 


0.02 


0.03 


Figure 5: *('0730} vs. symmetry-breaking mass m 3 at /3 = 1.0 on 24 3 x 20. The diamond shows the 
corresponding observable extrapolated to the chiral limit using staggered fermions 135] . 
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5.2 Meson Propagators 

Next we turn attention to correlations between mesonic operators ipY^ at different space- 
time points, with Y € { 75 , 73 ,1, ^ 7375 } all yielding spin-0 states, two with positive parity 
and two negative. We again focus on the recovery of U(2) symmetry as L s —> 00. A 
similar analysis for staggered fermions in the continuum limit of non-compact QED 3 was 
presented in pKJ]- We begin by listing two important identities of the 2+1+ld fermion 
propagator S(rrii ; x, s ; y , s') = (T(:r, s)^(y, s')), which follow directly from (1571) and have 
been checked numerically on small lattices: 

7 5 <S'(m/i, m 3 ,777.5; x, s; y, s')7s = S\m h , m 3 , -m 5 ; y, s'; x, s); ( 50 ) 

ld,S[rn h , m 3 , m 5 ; x, s; y, s')-/ 3 = S f (m h , -m 3 , m 5 ; y, s'; x, s). ( 51 ) 

Here x denotes the coordinate in 2+1 d and s that in the 3rd direction, s = L s — s + 1, and 
the dagger acts on Dirac indices. In what follows a mass parameter m, is omitted from 
the argument list if it takes the value zero. The identities (15015ip are of course crucial 
for the efficient calculation of timeslice correlators with a minimal number of matrix 
inversions, although there will be instances where a further inversion with the sign of 
the mass term m 3 $ reversed is needed. The contrast with relations (1551) for overlap 
fermions, which away from the continuum limit only hold for m 3 = m 5 = 0 should be 
noted. As will be demonstrated, in the limit L s —y 00 two further approximate relations 
become valid: 


'j 3 S(m 3 ,x,s;y,s')'y 3 « -S’ t (m 3 ; y, s'; x, s) for s' = s; 

lf 3 S(m 3 ; x, s; y, s')^/ 3 ?« ^^ 3 ; y, s'; x, s) for s' = s, (52) 


and equivalent ones for 3 +>• 5. 

Now consider the pion correlator C 3 (x) = ( , 0(O)75'0'0(ic)75t/ ; (^))- Under the parity 
definition dTO]) this state has J p = 0 _ . Using the definition (14T1) and ignoring the overall 
minus sign from the Grassmann nature of T, T, we deduce the following relation in terms 
of 2 + 1 +Id propagators 


C 5 {x) 


tr 


P + S{rrii, x, L s ; 0, L s )P_'y 5 P + S(m i ; 0, L s ; x, L s )P_^ b 


PS (mi] x, 1 ; 0 , l)P + ^ b P-S(mi, 0 , 1 ; x, l)P +75 


P_5'(m i ; x, l’,Q,L a )P_'Y 5 P + S(m i ]Q,L a ]x, l)P +75 
P + S(mi, x, L s ; 0, l)P + 7 5 P_S'(m i ; 0,1; x, L S )P _75 


(53) 


It’s natural in this formalism to have projectors P± sandwiching the fermion propagators 
- they are the bridge, both formally and in the code, between the Ad and 3d worlds. 
Now specialise to the case of mass term Sh, and use the identity (1501) . along with 
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translational invariance in 3d, to arrive at 



tr 


S(m h ; 0, L s ] x, L s )P_S\m h ; 0, L s ; x, L S )P + 


+ S{rrih, 0,1; x, 1 )P + S f (m h ; 0,1; x, 1 )P_ 

+ S(m h ; 0,1; x, L s )P~S\m h ; 0,1; x, L S )P_ 

+ S(m h - 0, L s ; x, l)P+S ] {m h ; 0, x, 1)P+ 


= C' /! - + (x) + C h+ “ (x) + (x) + P h++ (x). 


(54) 

(55) 


This has the familiar look of a pion correlator (ie. positive definite terms of the form 
SS^) apart from the projectors, and is clearly calculable with two sources, one located 
at s = 1 and the other at s = L s . In the abbreviated form (1551) note that —I- and H— 

correlators involve fermion propagators linking a domain wall to itself, whereas for- 

and ++ the propagators run between the walls. 

Next consider the 0 + correlator C 3 (x) = ( , 0 (O) 73 ' 0 ' 0 (x) 73 ' 0 (x)). This produces an 
expression with a different 4 d spacetime structure, but which with the help of (1511) can 
be rearranged to coincide exactly with (154)) : 


C&x) 


tr 


P-S(m h ] x, 1; 0, L s )P_'y 3 P_S(m h ] 0,1; x, L S )P_ 73 


PS(m h ; x, 1; 0, l)P + 'y 3 P + S(m h ; 0, L s ; x, P)P_ 73 


P + S(m h ; x, L s ; 0 , L s )P^ 3 P_S{m h ] 0 , 1 ; x, l)P+ 7 3 
P + S{m h ; x, L s ] 0, l)P + 'f 3 P + S(m h ] 0, L s \ x, l)P+ 7 3 
C h — (x) + C h+ ~(x) + C h - + (x) + C h++ (x) = C 5 ft (x). 


(56) 


These two states are thus exactly degenerate if the hermitian mass term S'/,, is chosen. 
The same methods can be used for the correlators 


C'i(x) = (^( 0 )^( 0 )^(x)^(x)); 

C 35 (x) = (^(O)i 7375 ' 0 (O)^(x)i 737 5 ^(x)), (57) 


which are respectively 0 + and 0 to show 

Cf (x) = C%(x) = tr S(m h ; 0,1; x, L s )-f 3 P_S ] (m h ; 0,1; x,L s )P _73 

+ S(m h ; 0, L s ] x, L s )^ 3 P_S\m h -, 0, L s ; x, L S )P+ 73 
+ S'(m fe ; 0,1; x, l) 7 3 P + S' t (m/ l ; 0,1; x, 1 )P_ 7 3 
+ S(m h ; 0, L s ; x, l) 7 3 P + S t (m/ l ; 0, x, l)P+ 7 3 

= C' h —(x) - C h - + (x) - C ft+ -(x) + C^x). 


(58) 


Therefore these two states are also degenerate, but distinct from i/j^ 3 ip and '075'*/’ i n 
(154)) . Assuming the primitive correlators C h±± are all positive dehnite, then \C£\ = 


17 








|G 3 | > \Ci \ = |C 3 5 1. If the symmetry breaking in the direction of S 3 is spontaneous 
this is consistent with U(2)—)-U(l)<g)U(l) yielding two light pseudo-Goldstone modes 
interpolated by 0750 and 0730. Had we instead chosen (HIT) as the definition of parity, 
then the assignments of the two Goldstone states would be reversed so that 0750 is 
now 0+ and 07 3 0 CT. However, the overall picture of a Goldstone O 31 doublet and a 
non-Goldstone CH doublet remains unaltered. 

The picture is more interesting when the mass term S 3 is considered. The correlator 
Gf = G 3_+ + G 3+ ” + G 3 + G 3++ in analogy to (j54M55j) . but for Cl we find 



tr 


S(m 3 ] 0,1; x, L s )P_S Jf (-m 3 - 1 0,1; x, L S )P_ 


+ S(m 3 ] 0 , l;x, l)P + 5’ t (-m 3 ; 0 , l;x, 1)P_ 

+ % 3 ;0,i s ;x, L s )P„S\ -m 3 ]0,L s -,x,L s )P + 


+ S(m 3 ; 0, L s ; x, 1 )P+5' t (-m 3 ; 0, L s ; x, 1 )P+ 


(59) 


= C 3 ~~ (x) + C 3+ ~ (x) + G 3_+ (x) + C 3++ (x). (60) 


The expression (160]) is similar in form to (156J) but requires twice the number of matrix 
inversions to evaluate. However, use of the large-L s approximations (|52|i gives 

Cl(x) « -C 3 —(x) + C 3 + ~(x) + C 3 _+ (x) - C 3 ++ (x). (61) 


Similarly we find 

C 3 (x) w -G 3 (x) - C 3 ++ (x) - G 3 _+ (x) - G 3+ "(x); (62) 

G| 5 (x) = G 3 (x) + G 3 ++ (x) — G 3_+ (x) — G 3+- (x). (63) 

This suggests that the Goldstone modes are now interpolated by ^ (0 + ) and ^ 75 ^ 
( 0 _ ), consistent with the variation of S 3 under (CMS]), but that the U(2)-»U(1)<8)U(1) 
pattern will now only be recovered in the limit L s —» 00 . 

Finally, for mass term S 5 we have 


Cl (x) w G 5_+ (x) + G 5+ - (x) - C 5 — (x) - G 5++ (x); (64) 

G|(x) = G 5 _+ (x) + G 5+_ (x) + G 5 (x) + G 5 ++ (x); (65) 

G 3 (x) = —G 5 _+ (x) — G 5 +_ (x) + G 5 (x) + G 5 ++ (x); (66) 

G| 5 (x) w —G 5 _+ (x) — G 5 +_ (x) — G 5 (x) — G 5 ++ (x); (67) 


the apparent contradiction with the expected identification of -0-0 (0 + ) and 0730 (0 _ ) 
(using either the symmetries ([21101) or the change of variables (1131) ) as Goldstone inter¬ 
polators will be further discussed below. 

We now present numerical results for the primitive correlators G ±=t (x). It is clear 
the recovery of U(2) symmetry hinges on the validity of the approximation (152|) . As 
in Sec. 15.11 the first test uses the numerically most demanding system of /3 = 0.5 on 
24 3 x L s , with fermion mass m l = 0.01. For a single equilibrated configuration we 
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Figure 6: Ratio J^fto \C/C 3± | vs. L s for C = C + (top) and C = C (bottom) for (3 = 0.5 on 24 3 . 

The key denotes the mass term used. 


calculated C + ~ and C , and where relevant their negative mass counterparts C ±_ , 
using five randomly chosen sources with s = 0. Using C 3± ~(L S = 48) as a reference 
point, the ratio Ylt=o C(L S )/C 3± ~{4&) is plotted as a function of L s in Fig. [6] for mass 
terms rrihSh, along with the same quantity evaluated for |(7 3 | (labeled “3-” 

in the plot). It is apparent that U 3±_ has reached its large-L s limit by L s & 40, but that 
C h ±- 

converges rather more slowly. Most importantly, C 3+ is indistinguishable from 
C 3+ ~, whereas C 3 actually changes sign around L s !=s 25 and by L s m 40 is practically 
equal to — C 3 , consistent with (l52]h The correlators C 5± ~ do not, however, fit the 
same pattern; C 5+ ~ lies systematically above C 3+_ , while C 5 —> 0 as L s —> oo. 
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Figure 7: Primitive timeslice correlators C 3+ ( t ) (left) and C h+ (t) (right) for various L s with /3 = 1.0, 

m = 0.01 on 24 3 . 



t t 


Figure 8: Same as Fig. |7]but for C 3 (t ) (left) and |C 3 (t)| (right). 

This picture is confirmed by a study of 500 configurations at (3 = 1.0, rn t = 0.01 
on 24 3 . Only the components C ± ~, which are manifestly real, are calculated; the coun¬ 
terparts C T+ are identical, and the correlators C ±T requiring an inversion with a neg¬ 
ative fermion mass recover these properties as L s —)■ oo. To focus on systematic ef¬ 
fects the same ensemble was used for each L s - the resulting statistical errorbars show 
the expected evidence of autocorrelation between timeslices. Fig. |7] compares C 3+_ 
with C hA ; C 3H converges to its large-L s limit for L s ^ 12, whereas C hA converges 
somewhat more slowly, requiring L s g) 16. Crucially though, by L s = 24 the two are 
indistinguishable. Fig. 0 confirms that after a sign change occuring for L s « 8 then 
lim^-Kx, C 3 = — C 3 , and that both also coincide in magnitude with C h for 
L s > 24. We have thus demonstrated that the mass terms Sh via relations (156158(1 . and 
S '3 via ((601163(1 . must yield equivalent meson spectra as —>■ oo. 
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Figure 9: Same as Fig. [7] but for C 5+ (t) (left) and C 5 (t) (right). Also shown is the combination 

C 3+ ~{t) + C 3 — (t) for L s = 24. 


Finally Fig. |9] shows the primitive correlators C 5± ~(t) for various L s . This time 
the picture is different: lim^^oo C 5+ ~ = C 3+ ~ + C 3 , while limx, s ^. 00 C' 5 = 0. The 

relations (I64H67P then predict that all four meson states become degenerate in this limit, 
and in this sense U(2) symmetry is realised. Explicit breaking to U(l)<g>U(l) by a mass 
term is not accomplished, however. We can trace this back to the expressions 

for the would-be Goldstones interpolated by 0730 (1651) and 00 (1661) which are exact 
even for finite L s , and can only coincide in magnitude in the limit that either C 5 or 
C 5+ ~ vanishes. There is thus no opportunity for degeneracy breaking induced by S 5 to 
develop; the system responds to this constraint by C 5H becoming larger as L s —> 0 so 
that all four mesons become degenerate with the two would-be Goldstones observed with 
mass terms mhSh,m^S 3 . Whether this exceptional behaviour is restricted to a singular 
point on the U(2) manifold or whether it is a more general phenomenon must be the 
subject of further study. 

6 Discussion 

The main focus of this paper is the numerical investigation of domain wall fermions for 
problems involving fermions in reducible spinor representations in 2+1 d. The results of 
Sec. 03 suggest it is highly plausible that the correct global U(2) symmetries are recovered 
in the limit L s -A 00 , independently of the continuum limit, just as for QCD [32]. How 
quickly the limit is approached remains, of course, a dynamical issue. However, the 
enhanced global symmetries of the problem admit new, helpful features not present 
in QCD, namely the ability to introduce a U(2)-breaking mass term in any of three 
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independent directions all yielding equivalent results as L s —> oo, aiding the robust 
determination of the L s required in practice. It is also significant that the dominant 
artifact A h(L s ), related to the “residual mass” in QCD simulations, cancels from the 
physical observables explored here when the twisted form im 3 , 5 ^ 73 , 5^ is used. 

The only exception to the general good news is the failure of the meson spectrum 
calculated with im^ 75 ^ to manifest the expected U(2)—)-U(l)®U(l) pattern, essentially 
because the formalism results in distinct expressions for the Goldstone correlators in this 
(possibly singular) case. The system responds by forcing their difference to vanish as 
L s — > 00 , so that no symmetry breakdown occurs. The general conclusion, however, 
must be that the method looks extremely promising, and optimal if symmetry breaking 
is implemented with a mass term im 3 7 / 73 -0. The next step is to explore a fully dynamical 
implementation of domain wall fermions, and then test U( 2 ) restoration by computing 
the separate components of the axial Ward identities. In this regard it is worth noting 
that the identity (150]) implies that the fermion determinant is real for mass terms Sh and 
S 3 , and hence there is no Sign Problem obstruction to running hybrid Monte Carlo with 
even Nf. For mass term S 5 the identity (l5Tj) does not suffice to prove detid D " real, but 
is completely analogous to the corresponding identity for domain wall fermions in 4 d. 

An interesting distinction with QCD, or with more general 4 d theories, is the relation 
between GW/overlap and domain wall approaches. As stressed, the two approaches 
differ in that GW formulation presented here maintains strict equivalence between 73 
and 75 whereas the domain wall treats them very differently. We have seen in Sec. [3] that 
for GW fermions the full U(2) symmetry is only recovered in the continuum limit: this 
is seen either via the effect of the remnant symmetry rotations 8 GW on bilinears ([Hip]) , 
or from the failure of the fermion propagator identity (1351) to generalise to 7713,5 7 ^ 0. 
Perhaps the problem arises from the GW implementation of the twisted mass terms (134|) . 
which as noted introduces an 0 (a) correction which cannot be compensated by a simple 
held rescaling. On the face of it, then, GW/overlap fermions can only recover U(2) for 
a —)■ 0, whereas domain wall fermions apparently recover U(2) for L s —)■ 00 irrespective 
of whether a continuum limit is taken. One might speculate that there is a different set of 
GW relations distinct to (I19II20II21|) . in which 73 and 75 do not appear on an equivalent 
footing, but which is consistent with the large-L s limit of domain wall fermions even 
away from the continuum limit. It would also be very valuable to understand things 
from a more fundamental perspective by attempting to derive the overlap from the large- 
L s limit of domain wall fermions, as done for 4 d in mi- It will also be interesting to 
extend consideration to the case of non-vanishing chemical potential / 1 , where there has 
been some discussion over the correct implementation of the GW symmetries [48, 49] : 
happily, there are several interesting 2+1 d systems where simulations at /i 7 ^ 0 with 
orthodox Monte Carlo methods are feasible [snusmsa. 

The biggest goal, however, once dynamical fermions have been implemented, is to 
study the strongly coupled fixed points of the GN, Thirring and graphene [53] [54] systems 
to see if the critical exponents match those previously obtained with staggered fermions, 
and to check whether the expected dependence on Nf is seen (an exploratory study of 
the GN model with domain wall fermions is presented in [55] )• This would not only 
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provide an important and hitherto unexplored test of the applicability of the domain 
wall approach to lattice fermions, but more generally would represent an important step 
in the definition of interacting quantum field theories beyond weak coupling. 
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